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We investigate the binding energy in two classes of poly- 
tropic perfect fluids. A general-relativistic bound from below 
is derived in the case of a static compact body, having the 
same form as in the newtonian limit. It is shown that the 
positivity of the binding energy implies that the (properly 
defined) average speed of sound is smaller than the escape ve- 
locity. A necessary condition for the negative binding energy, 
stating that the maximal speed of sound is close to the escape 
velocity, is found in a class of fluids. 



I. INTRODUCTION 

There are two popular classes of equations of states of 
perfect fluids. One, that is commonly used by astrophysi- 
cists, takes a simple form when formulated in terms of 
the so-called baryonic number density. The other class, 
preferred by general relativists, uses the mass density. 
The two classes coincide in the newtonian limit but they 
can lead to markedly different predictions in the strong 
field regime (see, for instance, a comparative discussion 
of stationary accretion in [1] ) . ft is intention of this pa- 
per to rederive some known results and obtain new ones 
on the binding energy of perfect fluids satisfying poly- 
tropic equations of state of the two types. A possible 
new feature of the derivation is that we do not recourse 
to the laws of relativistic thermodynamics. We restrict 
our attention to static fluid bodies, but since we operate 
in those classes of equations of state for which staticity 
implies spherical symmetry [2], we end in investigating 
spherically symmetric objects. 

The order of this paper is as follows. Section 2 presents 
equations of selfgravitating fluids. In the third section we 
recall definitions of the two classes of relativistic equa- 
tions of state and show that in the (properly defined) 
Newtonian limit they coincide. Section 4 defines the 
gravitational binding energy and in Sees 5 its value is 
estimated in the polytropic subclass of the equations of 
state. In the newtonian limit the binding energy agrees 
with the standard newtonian expression. Section 6 is 
dedicated to the derivation of certain inequalities onto 
the gravitational potential energy, ft is remarkable that 
these inequalities are exact, and one of them has the same 
form both in the general-relativistic and in the newtonian 
cases. Section 7 presents the main results. It is shown 
that the binding energy of a compact spherical body is 
bounded from below by the standard newtonian formula 
-P/(r - 1) + 3m^ B /(5R). The kinetic energy of the gas 
is represented by the total pressure P and the sufficient 
condition for having a positive value of the binding en- 



ergy is that P is smaller than the gravitational potential 
type term 3m^ B /(5R). That can be translated into a 
bound from above onto a quantity that represents a kind 
of an average speed of sound - that should not exceed the 
escape velocity (or even the first cosmic velocity). A nec- 
essary condition (a lower bound on the maximal speed of 
sound) for the negative binding energy is formulated in 
the so-called p — p class of fluids. 



II. EQUATIONS FOR SELFGRAVITATING 
FLUIDS. 



Throughout this paper the Newtonian constant G and 
the velocity of light c are put equal to 1. Let the self- 
gravitating system be spherically symmetric. Then the 
whole dynamics is carried by a material field character- 
ized by an energy- momentum tensor Tp„. fn the case of 

Pg^v, where the 
— I. 



the perfect fluid = (p + p)U^U v 
coordinate velocity is timelike, U^U^ 

Assuming spherical symmetry, one can deal with the 
metric 



ds z 



-N 2 dt 2 



adr 



R 2 dn 2 , 



(1) 



where N, R and a are unknown metric functions. In 
sections 2, 3 and in the first part of Sec. 4 we work in 
comoving coordinates. Our aim there is to justify a par- 
ticular expression of the binding energy, and that requires 
some information about the evolution that led to the for- 
mation of a bound system. Comoving coordinates are 
convenient in describing evolving fluids. Starting from 
section 4 we use standard (polar gauge) coordinates. 

The needed material quantities are the energy density 
p = —Tq and the isotropic pressure, p = T£ = T[j = Tt. 
Define pR = 2d r R/y / a (this is so-called mean curvature 
of a centered sphere [4]) and U = dtR/N. One can show, 
using the constraint equations [3] and the properties of 
comoving coordinates, that 



m(R) 

/'/'' = 2\/l — 2 — jj— + U 2 , 



(2) 



where m(i?) is the quasilocal mass specified below. The 
trace of the extrinsic curvature 



trK = K] 



jjd In ^det(g. 



can be found from the momentum constraint [3], 

M ur2 ) 



trK = 2- 



^/ER 3 p 



(3) 



(4) 



1 



The continuity equations ViTj = reduce to two equa- 
tions 



Nd r p + d r N(p + p) = 0, 



d p = -NtrK(p + p). 



(5) 



(6) 



The Einstein evolution equations reduce to the single 
equation 

d U = -N^ + d -I^-^Nm. (7) 
R z a 

The quasilocal mass (" mass function" ) 

>(r) = I dV^-p (8) 

JV(r) 1 



contained within a coordinate sphere r' < r satisfies the 
evolution law 



d m(r) = -4irNR 2 Up. 



(9) 



The mass m(r) (which happens to be equal to the vol- 
ume representation of the Hawking mass [4]) can be ex- 
pressed also as 



m(R(r)) = 4-7T 



R(r) 



drr p; 



(10) 



the integrand is now written in terms of the areal ra- 
dius. The total mass is conserved; the asymptotic mass 
to(oo) is constant on the whole foliation that develops 
from a given initial hypersurface. There exists another 
mass measure, known as the rest mass [5], 



M(R) 



dVp = 8tt 



V{R) 



f 

Jo 



dr-p. 
P 



(11) 



The rest mass M(oo) is not conserved - it might change 
from a slice to a slice (that is with time). In what follows 
we will write M and m instead of M(co) and m(oo) 

The baryonic mass (the mass density n is defined in 
the next section) contained within a coordinate sphere r, 



m bar (r) = / dVn 

'V(r) 



r rR(r) 
/ dVn = 8TT I 

Jv(r) JO 



dr-n, 
P 



(12) 



is conserved: domi, ar (r) — 0. That follows easily from 
equations (16), (6) and (3). Yet another mass charac- 
teristic of a static body comprised within the coordinate 
sphere r is its Trautman-Bondi mass ttitb, equal to the 
quasilocal mass function m(r). This mass is equal to the 
asymptotic mass minus that fraction of the fluid that has 
been dispersed away to infinity. Therefore ttltb < rn. 



III. EQUATIONS OF STATE 

The (p — n henceforth) equation of state used by as- 
trophysicists is 



p = p{n) 



(13) 



while the community of general relativists prefers (p — p 
thereafter) 



p = p{p). 



(14) 



Here n is the baryonic mass density that can be defined 
(up to a normalization factor) as an integrability factor 
that ensures the conservation of the baryonic current 



(15) 



One can easily show [3] that the conservation laws 
V M X^ = imply, under conditions assumed above, the 
relation 



n = hq cxp 




1 



■p(s) 



(16) 



Here n and p are constants that can be used in order to 
suitably normalize n. Thus, given p and an equation of 
state p = p(p) satisfying V^T^ 1 = 0, one finds n automat- 
ically satisfying the continuity equation V At (n[/ M ) = 0. 

And conversely, given conserved baryonic density and 
an equation of state p = p(n) one can find p, 



ds 



p{s) 



n, 



(17) 



such that V M T^ 



0. The formulae are particularly sim- 
ple in the case of polytropic equations of state. Assuming 



p = Kp l 



with T > 1 one has 



n - 



n - 



(1 + f) 5 



(18) 



(19) 



(a is the speed of sound). If one assumes no = 1 then, 
defining the newtonian limit to be the case a 2 <C 1, it 
follows that in the newtonian limit the baryonic density 
and the mass density are very close, n « p. 
Assuming the alternative state equation, 



p = Kn r , 



one obtains from Eq. (17) 



(20) 



(21) 



if the newtonian limit is understood as the stronger than 
before condition p/(T — 1) <C 1, then d p n s=a 1 and again 



2 



n w p. Assuming in Eq. (19) that no — > and simultane- 
ously keeping the pressure p fixed (that is KtIq = const) , 
one obtains the Harrison et al. equation of state [6] 



P = 



P 



r- 1 



(22) 



In the following part of this work we shall deal with 
equations of state given by (18) and (20). The adiabatic 
index T is assumed to be strictly bigger than 1. While 
the above relations are well known, it is interesting that 
the above derivation can be done solely within classical 
general relativity, without appealing to thermodynamics. 



IV. BINDING ENERGY 

The binding energy E of a static system S can be de- 
fined as the difference between the initial mass that is 
needeed in order to create S and the final Trautman- 
Bondi mass. The most natural definition would be to 
put E = m — niTB: E is the difference between the 
asymptotic (ADM) mass m and the mass wtb of the fi- 
nal product. Traditionally one defines the binding energy 
somewhat differently, using the concept of the baryonic 
mass. We will adhere to the standard approach, defining 
another measure Ebind of the binding energy; the reader 
may check later that E > Ebind- 

In the case of initial mass distributions that are at 
rest and so dispersed that the curvature of the Cauchy 
hypersurface is practically negligible, the two masses m 
and M can be practically equal. If we assume that matter 
consists exclusively of baryons, then (again, under the 
condition that the fluid is so diluted that its pressure 
is negligible) the baryonic mass mf, ar (r) and the mass 
function m(r) can be regarded as being equal initially. 
Since usually only a fraction of the gas (say, comprised 
within the coordinate sphere r) becomes solidified in a 
form of the static object S, and it is the baryonic mass 
mbar(r) that remains constant during a collapse, it is 
justifiable to define the binding energy E as [7] 



Ebind = m bar (r) - m T B- 



(23) 



Let us point out that the baryonic mass can be expressed 
- using formulae of Sec. 3 - in terms of quantities de- 
scribing the static system S. From now on we will employ 
standard (Landau-Lifschitz [8] ) coordinates with the ra- 
dial coordinate being always the areal radius. This im- 
plies a = 1/(1 — 2m W ) , where the mass function is given 
by Eq. (10). All volume integrals will be performed over 
the volume V of S. In order to simplify notation, the 
total rest mass M(r) of S will be denoted by Ms and the 
binding baryonic mass TOfc ar (r) will be written simply as 
nib ar . Notice that totb is this quantity that enters the 
newtonian formula describing the late-time interaction 



of S with a distant (but located at a finite distance) test 
body. 

Let us define total pressure in S by 



P = [ dVp = 8n [ dr-p. 
Jv Jo P 



(24) 



It is convenient to rewrite formula (23) in two different 
forms, 



Ebind = M S 



p 



r - 1 



m TB - 



in the case of the p — n equation of state and 



Ebind = / dV- — — ; m TB . 

v (l + £)r=r 



(25) 



(26) 



p — p equation of state. The equivalence of (26) and (25) 
with (23) follows from (21) and (19), respectively. 

In the newtonian limit both expressions yield the same 
quantity 



V. ESTIMATES OF THE BINDING ENERGY 

One can invoke to a workable version of the virial theo- 
rem, that exists in static spherically symmetric systems, 
in order to eliminate P from (25) ( [9] - [11]). This is 
obtained from the Oppenhcimcr-Volkov equation 



{p + P )(^ + ^pR)=-d R p(l-^}, (28) 



2m{Ry 



Multiplying (28) by 4irR 3 / \J 1 - and integrating 

with respect R, one gets (integrating by parts, in order 
to eliminate the grad-term dnp) 

-U g + Att / dVpr 2 (p + 2p). (29) 
Jv 



Here Ug is the gravitational self-interaction energy, 



U G 



L 



dV 



pm(r) 



(30) 



The total pressure P = J v dVp is not smaller than P, as 
one can see rewriting the first line of (29) 



P = P+ [ dV 
Jv 



2m(r)p 



(31) 



The mass difference Ms—niTB reads (this follows directly 
from definitions of the two masses) 



3 



M s - m TB = -U G 



Define 



dV- 



2pm(r) 



2m(r) 



(32) 



■to—Oo + WTify 



dV- 



2pm(ff 



r 2 I + 



2m(r) 



2 - 



(33) 

Employing the above information, after some calcula- 
tions one arrives at 



M s 



r - 1 



- rriTB = 



M s 



— - mTB -r 



hi 



dV 



2m(r)p 



3r-4 



U G 



3(r 

Thence, 



3(T-1) " (T-l)jv 
4 ^y J dVpr 2 (p + 2p) . 



2m(r)p 



p sr-4 - 



r - 1 



3(r-i) 



(34) 



(35) 



Notice, that Eu-nd — Ms — fAj- ~ m TB for the p — n 
equation of state. Therefore we infer that the binding 
energy may become negative for values of the polytropic 
index T < 4/3. In the newtonian limit one obtains 



Ebind - ~ 3(f _ 1) ^ G " 



(36) 



In the newtonian limit the binding energy becomes 
strictly negative only when T < 4/3, while in the regime 
of strong gravitational fields this may happen even for 
T = e + 4/3, for small but positive e. The general rela- 
tivistic effects may decrease the binding energy (in agree- 
ment with [7]). Let us remark here, that the result in [7] 
corresponding to Eq. (36) has a different coefficient (for 
reasons not quite clear to us): (3r — 4)/(5T — 6) instead 

of (3r-4)/(3(r-i)). 

Notice that always E^nd < Ms — vhtb- In the nonrel- 
ativistic case and putting T = 5/3 one obtains a stronger 
result, Ebind = —Uq/2 w (M — totb)/2; the quantity 
Ms — tub, that can be regarded as the gravitational po- 
tential energy in the case of static systems (or perhaps 
even momentarily static systems), bounds the binding 
energy and it is of the order of the binding energy. The 
second part of the last statement is not valid for strongly 
curved geometries, even in the case T — 5/3. 

Assuming the polytropic p — p equation of state, we 
can express the binding energy by formula (26). It is 



easy to show that y J > 1 — ffpirq ■ Tlicrctoro 
one arrives at 



Ebind >M S - 



P 



r- i 



m TB - 



(37) 



It is clear that now the use of the virial theorem (which 
bounds from above M s — p^y — m TB ) cannot give an 
estimate from above on the binding energy. There ex- 
ists only a trivial bound Ebind < — Ug (this is because 
Ms > rribar) - the binding energy is bounded by the 
absolute value of the gravitational potential energy. One 
can see, comparing formulae (34) and (37), that the p— p 
equation of state is more favourable for having nonneg- 
ative binding energy that the p — n case. Nevertheless, 
in the newtonian limit one again recovers estimate (36). 
We will show in the last section that by exploiting the 
quantity M — vtitb one can find another useful bound 
from below. 



VI. BOUNDING POTENTIAL ENERGY 

In this section one might partly relax former con- 
ditions, and to deal with a compact (not necessarily 
fluid) body on the momentarily stationary spatial 3- 
dimensional geometry. The assumptions i) of staticity, 
and ii) that matter is a barotropic fluid, p = Kp T ( p \ 
become relevant in a refined version of the forthcoming 
result. To be specific, the second and third terms on the 
right hand side of Eq. (39) would be absent without the 
two last conditions. Let us point that the compactness 
of the fluid body imposes a bound from below on the 
barotropic index ( [12] and [13]). 

It is easy to show the following result, that estimates 
the difference Ms — i«tb • 

Theorem. Assume that matter has compact support 
enclosed within a sphere R, that no minimal surface does 
exist on the Cauchy slice and p > 0. Then 



-Ug = M s - mTB > 



rn 



TB 



R 



(l + V 1 -^) 



> 



rn 



TB 



2R ' 



(38) 



If in addition the configuration is static, then 



-U G = M s - rriTB > 




R . , / 2m T B N \ . 

arcsin(W— — ) > 



(39) 
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Proof of the theorem. Straightforward calculation al- 
lows one to check that 

1 f R ( - d 1 \ 

Ms = 77 / dr\ -2^?—^r-2m{r) + = 

2 Jo \ dr w L 2m(r) y 



il l 



2m(r) 



The integrand of the last integral is bounded from below 
by 2. Therefore M s > R(l - - ^ s -); thus 



M s - m TB > R(l 2m I B ) - m TB 



R 



l TB 



R 



> 



rn 



TB 



2R 



(41) 



We can get only thus far in the case of momentarily sta- 
tionary data. 

If our configuration is static, then the Oppenhcimer- 
Volkov equation implies dup = a~ 2 dnp < 0. This in 
turn yields d r (m(r)/r 3 ) < (that result was originally 
derived by Buchdahl [14]). 

Returning to Eq. (40), notice that 



1 



ft 



dr 



1 



R + 



R dr 



1 



2m 1 



2m(r) 



o 



(i+ymf)' 



(42) 



Since m(r)/r 3 is decreasing, we have m(r)/r 3 > 
rriTB/R?- Therefore 



f n dr m 2 
Jo ^^1-^(1 + ^/1-^) 

I 



> 



2m, 
~R 6 



2 r R 

TB 



1 + 



2m,TBr 2 



, (43) 



The second integral term in Eq. (43) can be explicitly 
evaluated, with the result 



F = 



-mrB 




R . 2mTB 

arcsm 1 



2m TB 



R 



(44) 



Combining (40), (43) and (44), one gets the first 
of the sought inequalities of (39). The second 



inequality there follows from the trivial estimate 
4 J R dr ; . rA : > R 5 /5. That ac- 



R3 



R3 



complishes the proof of the theorem. 

Let us again stress here that while in (38) one needs 
only that initial data are momentarily static, the result 
(4Qfhth the coefficient 3/5 in (39) requires staticity (with- 
out that we would have only the coefficient 1/2). The 
two results are exact. The first inequality is saturated 
by the spherical shell of matter (see an example stud- 
ied in [15]) while the inequalities of (39) are saturated 
both in the general-relativistic and the newtonian cases, 
by a constant density star. In the case of fluid bodies 
satisfying the Buchdahl limit R = QrriTB/^ [14] one ob- 
tains M s > (-9/16 + 81 arcsin(2V2/3)/(32V2)) m TB « 
1.65raTs; that again becomes an equality for constant 
density stars. Constant density stars are the most eco- 
nomic configurations in the sense that they minimize the 
ratio Ms /mTB - and, consequently, the gravitational po- 
tential energy —U- amongst the subset of static solu- 
tions. 



VII. CONCLUSIONS 

Before pursuing further, one should be aware that the 
assumption of compactness of a fluid body made in Sec. 6 
imposes some restrictions on the polytropic index. In the 
case of p — p equations of state the condition is roughly 
T > 6/5, but for the p — n case the situation is unclear - 
the application of existing results concerning the general 
barotropic models ( [9], [12] and [13]) requires a separate 
investigation. We assume that the polytropic index is 
just right to agree with the compactness. 

Collecting together the results of Sees 5 and 6, we ar- 
rive at the inequality, valid for both p — n and p — p 
models, 



Eund > M s 



rriTB > 

,2 



r- 1 



+ F + 



rn 



TB 



R 



(45) 



A weaker but more elegant than (45) is the inequality 

(46) 



TP -> P 4- 3rrLl 

r^bind ^_ r— 1 



5ft ' 

which retains the same form in the newtonian limit. The 
sufficient condition that the binding energy be positive is 

^ 2 



P r-, rn'TR 
< F+ 1B 



r- 1 



R 



(47) 



or, more simply, 



3m TB P 



5R 



r- 1 
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both for p~n and p— p equations of state. An alternative 
form of a sufficient condition, 

P Mj 

r - 1 < ~2L 

can be obtained due to another inequality Ms — rriTB > 
Mg/(2L) (where L is the geodesic radius of a com- 
pact body), proven in [15]. There exists a simple in- 
terpretation of this condition in terms of the speed of 
sound. Namely, define an "average speed of sound a " 
by P = a 2 M s /T. The quantity a is bounded from be- 
low and above by the minimal a m in and maximal a ma x 
speeds of sound, respectively, within the compact body. 
The formerly formulated sufficient condition implies now 



r,2 



<r(r-i) 



3m; 



TB 



5RM S 



But Ms > rriTB', therefore 



*<r<r-i)^ 



(48) 



(49) 



as the consequence of the positivity of the binding energy. 
In the newtonian limit Mg « tutb and (49) becomes just 
another sufficiency condition for Eund > 0. This speed 
of sound should not be bigger than (roughly) the es- 
cape velocity (or, strictly saying, the first cosmic velocity 
ttltb/R), if r < 5/3, and it decreases with the decrease of 

3m 2 

T. And conversely, the inequality (a) 2 > T (T — 1) ^Jf s 
is necessary for Ms — jAr — rriTB being negative, and 
thus necessary for the binding energy being negative. 

Another form of the necessary condition can be found 
for p — p equations of state. It follows from (26) that 



Ebind 



-L 



dVn — rriTB > 



M s 



(1 + %-)' 



m TB . (50) 



The use of the estimate (39) of the Theorem allows us to 
say that Eund can be negative only if 



/ 



r-i 



n > r 

"•max — i 

-r 



1 + F + 



rriTB 



V 



R 



(i + V / i 7 ¥)' 



(51) 



or (in a simpler and weaker form) 
a 2 > r 



'IjutbV- 1 
5R J 



(52) 



Since x = < 4/15, one has (1 + x) S > 1+5 (±§) 

Therefore Eq. (52) can be written as 



5-1 



>|r,r 



\19 



2-r 



m TB 



(53) 



This clearly demonstrates, that the negative binding en- 
ergy requires the maximal speed of sound to be close to 
the escape velocity. For compact objects satisfying the 
Buchdahl limit we get from (51), putting V = 5/3, the 
necessary condition a 2 ^^ > 0.66. Let us remark that 
(52) and (53) yield much worse estimates a? max > 0.28 
and o 2 naa . > 0.27, respectively. In such a case the up- 
per limit for the maximal allowed kinetic energy (with 
Ebind > 0) exceeds (T — 1) x 0.65totb (see the end of 
the preceding section); the kinetic energy can constitute 
a significant fraction of the mass mj b ■ 

By reversing the argument, one obtains from the con- 
verse of (51) - (53) a set of sufficiency conditions for the 
positivity of Eund- The strongest statement is that if the 
maximal speed of sound is smaller than the right hand 
side of (51), then Eund > 0. 

As is well known, the negative sign of the binding en- 
ergy is correlated with the instability of static fluids [7]. 
Thus inequalities (48) and (51 - 53) agree well with the 
intuitive notion that hot systems (i. e., with high speed 
of sound) arc unstable. It is also remarkable that the 
upper limit increases monotonically with the increase of 

r. 

We conjecture that inequalities (45 - 53) hold true for 
systems with the general barotropic equation p = Kp v ^ 
(under a suitable definition of the constants). A plausi- 
ble form of a related condition in the case of stationary 

nonsphcrical fluid bodies can be —zr[ < m TB I \J~^ n 
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